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1 Introduction 



The study of the multiplicity distribution of the produced hadrons along with the analysis of 
the correlations among them stands in the frontier of investigations in the area of multiparticle 
dynamics. The multiplicity distribution plays a fundamental role in extracting first information 
on the underlying particle production mechanism, while the correlations give details of the dy- 
namics. Whereas the full-multiplicity distribution is a global characteristic and is influenced by 
conservation laws, the multiplicity distributions in restricted phase space domains contributing 
to the correlations are local characteristics and have an advantage of being much less affected 
by global conservation laws. 

In the last decade, study of multiplicity distributions in limited regions (bins) has attracted 
high interest in view of search for local dynamical fluctuations of an underlying self-similar 
(fractal) structure, the so-called intermittency phenomenon. For review, see Refs. |E0|. This 
phenomenon is seen in various reactions; however, many questions about intermittency and, in 
particular, its origin are still open and further investigations are needed [fj. [| . 

Studying local fluctuations, one must remember that the fluctuation of a given number of 
particles, q, is contributed by genuine lower-order, p < q, correlations. To extract signals of 
these p-order correlations, one needs to use advanced statistical techniques such as normalized 
factorial cumulant moments (cumulants) @-^[. However, this method requires measurements to 
be of high statistics (and of high precision), lack of which leads to a smearing out of high-order 
correlations @,^|. 

A real opportunity to extract genuine multiparticle correlations came with vast amount of 
multihadronic events collected now at LEP1. The statistics available allows one to decrease 
significantly the measurement errors and to reveal relatively small effects. The first reports 
representing the studies of genuine correlations in e + e~ annihilation have just appeared HQ- 
DELPHI has analysed correlations in one- and two-dimensional angular bins of jet cones, 
while OPAL [[8] has performed its study in three dimensional phase space using conventional 
kinematic variables such as rapidity, transverse momentum and azimuthal angle. DELPHI has 
shown an existence of the correlations up to third order, and OPAL with increased statistics 
has, for the first time in e + e~ annihilation, calculated multidimensional cumulants and has 
established sensitive genuine correlations even at fifth order. Note that multidimensional anal- 
ysis carried out in hadronic interactions shows that the cumulants are consistent with zero at 
q > 3 there, while in heavy-ion collisions non-zero cumulants have been observed at second 
order only 

The genuine correlations measured in e + e~ annihilations are found to exceed considerably |7| 
the QCD analytic description and to indicate || significant deviations from Monte Carlo (MC) 
predictions. A smallness of genuine correlations predicted by perturbative QCD is seen also in 
cumulant-to-factorial moment studies, even when higher orders of the analytical approximation 

1 Earlier, using correlation (strip) integrals to reduce statistical errors, multiparticle genuine correlations have 
been searched for in hadron-hadron and lepton-hadron |llj interactions. It is worth to note that except for 

problems arising due to various possibilities in defining of a proper topology of particles and a distance between 
them (see e.g. ||), a rather important difficulty in interpreting results could come from a translation invariance 
breaking of the many-particle distributions This will lead to different results on moments/cumulants 

according to a variable used because different variables are sensitive to different hadroproduction mechanisms 
|T^,|l3|. For example, in the high-order genuine correlations obtained in Ref. ]9|,[ll[] the study is performed in 
four-momentum difference squared, Q 2 , dependent on Bose-Einstein correlations, whereas in the pseudorapidity 
analysis no correlations higher than two-particle ones were found, since (pseudo)rapidity seems jl3j to be 
more "natural" variable to search for jet formation than e.g. for Bose-Einstein correlations. 



are used fn]>0]- These findings along with the deviations obtained by L3 in its recent detailed 



MC analysis of local angular fluctuations [|T^| show particle bunching in small bins to be a 
sensitive tool to find differeneces between parton distributions treated by QCD vs. hadron 
ones detected in experiments. The study of bunchings seems to be critical to a choice of a more 
convenient basis for a suitable approach of multiple hadroproduction. 

In this paper we compare the intermittency and correlation results from OPAL with pre- 
dictions of various parametrizations, or regularities. We consider the negative binomial distri- 
bution, its modified and generalized versions, the log-normal distribution, and the pure-birth 
stochastic production mechanism. For the first time we examine these parametrizations with 
high-order genuine correlations. The incorporation of the multiplicity distribution in the study 
of correlations provides more advanced information by using various approximations and mod- 
els. In particular, this study gives more understanding about the structure of multiparticle 
correlations, e.g. their relation to two-particle correlations. 

All the above listed parametrizations are well-known in multiparticle high-energy physics 
for a long time and are used to describe the shape of the multiplicity distribution, either in 
full phase space or in its bins. Essentially all these parametrizations are sort of branching 
models and show an intermittent behavior. This feature becomes particularly important in 
hadroproduction studies in e + e~ annihilations where parton showers play a significant role and 
the hypothesis of local parton-hadron duality is applied to the hadronization process M, O . 



2 Normalized factorial moments and cumulants 
from various parametrizations 

There is a variety of models which describe particle production as a branching process, see 
e.g. |r7J . The main prediction of these models is a suitable parametrization for the multiplicity 



distribution. Further, more details of the underlying dynamics come from investigations of the 
factorial moments and cumulants of the multiplicity distribution given ||. 

One of the most popular parametrization used to describe the data for full and limited phase- 



space (basically, rapidity) bins is the negative binomial (NB) distribution, see Refs. | IT IS 
for review on the subject and its historical development. 

The NB regularity depends on two parameters: (n) is the average multiplicity and 1/k is 
the so-called aggregation coefficient which influences the shape of the distribution. For the NB 
distribution, the normalized factorial moments and cumulants are governed by the parameter 
1/k and can be derived RPUf from the following formulae, 



F q = F q ^[l + q -l) (1) 



and 

K q = {q-l)\k l -\ (2) 

respectively. 

Fractal properties of the NB distribution have been studied elsewhere P,|21|-p4| . To demon- 



strate the fractality of the NB distribution, it has been proposed to use its generalization |22 
or to study individual sources obtained by consideration of different topologies, e.g. number 
of jets |24]]. A good fit by the NB parametrization should yield k independent of g-order, but 
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the limited statistics may also stabilize the g-dependence and hide the (unknown) real distribu- 
tion |j2l |. The bin-size dependence of the k parameter, expected from the bin-size dependence 
of the factorial moments and cumulants, Eqs. ([I]) and (0), is a consequence of the unstable 
nature of the NB distribution, i.e. the convolution of the distributions in two neighbouring bins 
does not give the same type of distribution p|. In the other words, the bin-size dependence 



of k reflects the fact that distribution in one bin depends on that in the other bin . 

In e + e~ annihilations at the Z° peak [[Z7|-|3"T] and lower energies [32] , the NB parametrization 
was found to fail consistently in describing the multiplicity distribution, either in rapidity bins 
or for the full phase space. Using e + e~ results on intermittency PSJ, it was also shown j33| that 



this parametrization does not reproduce the large fluctuation patterns, while it is appropriate 
for phase space bins in which the fluctuations are sufficiently small. Such an effect has been 
observed in (pseudo)rapidity studies in other types of collisions fl20|, |35H38| toofl. The likely 
reason is the above noted instability of the NB distribution. Indeed, the NB law provides an 
acceptable description in the central regions of the multiplicity distribution, away from the 
tails [17|, ^] . In this region, the distribution measured is mostly flat and, therefore, less 
sensitive to instability effects. 

Another reason is that the NB regularity underestimates the high multiplicity tail |p4|, |36|, 
pH^ . -EU, ¥1] which gives the main contribution to the fluctuations and which is influenced by 
instabilities. For high multiplicities the NB distribution transforms to the stable T-distribution 



25| . This type of distribution was found |27| to be the most adequate to describe the multiplicity 
distribution of the OPAL data. This is in contradiction to the results obtained in nuclear 
collisions |38| , where the T-distribution was found to be significantly inconsistent to reproduce 



the measurements: it underestimates the low-multiplicity parts of the experimental multiplicity 
distributions in different rapidity bins, while overestimates the high-multiplicity tails. Again 
in contrast to e + e~ data, the NB regularity is found to be the best one to describe small 
fluctuations in the multiplicity distribution in nuclear data, and large fluctuations are well 
reproduced by two-particle correlations |3l||¥U|-pE3]]p|. 

Another popular choice for the parametrization of the multiplicity distribution is the log- 
normal (LN) or Gaussian distribution ||44l,|45|l . This type of regularity of final state distribu- 



tion can be obtained by assuming a scale invariant stochastic branching process to be the basis 
of the multiparticle production mechanism. The LN distribution is defined by two parameters, 
the average and dispersion. To describe the data a third parameter has been introduced [EE 



to take into account an asymmetry in the shape of the full phase space distribution measured. 

In e + e~ annihilation this model has been found to describe successfully the data for the full 
rapidity window |27],[30l . For restricted bins the best agreement between the LN parametrization 



and the data has been obtained for very small bins in the central rapidity region, while for the 
intermediate size bins the deviation observed is assumed to arise from perturbative (multi-jet) 
effects 



30 



Compared to the NB, the LN regularity is found to give much better description 
which leads to understanding the multiparticle production process as a scale invariant stochastic 
branching process. 

The fluctuations have been studied in a model of this type in Ref. pi, within the so-called 



'a-model" and it was found that in this particular case the normalized factorial moments obey 



Comparing the formation of dense groups of particles in e + e~ and in hadronic (or nuclear) interactions, 
one finds noticeable difference of the fluctuations structure being isotropic (self-similar) in the former case and 
anisotropic (self-afhnc) in the latter one |}9| . This is expected to reflect different dynamics of the hadroproduc- 
tion process in these two types of collisions. 
3 See preceding footnote. 
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the recurrence relation, 



l„F, = Mi_i) |nF2 . (3 ) 

Strictly speaking, this formula connects the standard normalized LN moments C q = (n q )/(n) q 
rather than the factorial moments F q . The difference between this two types of moments 
becomes negligible at large multiplicities which is not the case for small bins. Thus, the re- 
sults based on the Eq. (|3|) with F 2 defined by data can deviate from the true LN regularity 



predictions, particularly at small bins |p4 | 



Recently, the modified negative binomial (MNB) regularity has been introduced to 



correct deviation between the NB parametrization predictions and the e + e and pp data E7 



One finds f 4S 1 for the normalized factorial cumulants of the MNB distribution, 



Here, r = A + (n) / k and the superscript minus indicates that this law is applied for negatively 
charged particles. The MNB regularity reduces to the NB one if A = 0, cf. Eq. (0). 

The MNB parametrization has been found to give an accurate description of the full phase- 
space multiplicity distributions in e + e~ annihilation measured from a few GeV up to LEP2 



energies |2q,f44f49|-|52|l and in lepton-nucleon scattering data in the wide energy range [p3l ■ A 
similar energy dependence of the parameter k has been obtained in these two types of collisions. 
Recently, a simple extension of the MNB law has been found to describe the charged particle 
multiplicity distributions in symmetric (pseudo)rapidity bins as well [Q . 



In contrast to the NB, the MNB parametrization is shown to reproduce fairly well the 
factorial moments and the cumulants of the full like-sign, e.g. negatively charged particle 
phase space from e + e~ data at the energies ranging from 14 to 91.2 GeV The fits of 

the multiplicity distributions, the moments and cumulants give rise tofc>0,A<0Q and 
0<r<|A|<lQ, that, according to Eq. ([|), leads to the cumulants being negative at even 
values of order q and positive at odd ones. This fact has been utilized to explain the oscillations 
of the ratio of the cumulants to factorial moments as a function of q . 

To obtain the quantities for all charged particles one uses the fact that the number of charged 
particles produced in e + e~ collisions is twice that of negative ones. Then, the normalized 
factorial moments and cumulants are given [48[] by 



The stochastic nature of the NB and MNB laws allows to generalize them by introducing 
a stochastic equation of the pure birth (PB) process with immigration |58|,|59[. Then, the 
NB and the MNB distributions can be derived from this birth process under the appropriate 
initial conditions, namely, the birth process with no particles in the initial stage leads to the 



4 Negative values of the parameter A are interpreted as the probability —A of intermediate neutral cluster 
to decay into the charged or neutral hadron pairs 49, 5(J. Positive A values are also acceptable Jt7| ) |5q] ) but 



in this case k is the maximum number of sources at some initial stage of the cascading [56] changing from one 
event to another, in contrast with the fixed k value corresponding to the case of A < [49 1. 

5 Nevertheless, at LEP1.5 energy, -y/s ~ 133 GeV, an inverse inequality |A| < r ~ 0.91 can be found pip . 
An increase of the parameter r with increasing energy is allowed but this complicates the particle-production 
scenario leading to the MNB parametrization p7|. 
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NB law, while the MNB distribution is resulted from the birth process with the initial binomial 
distribution |H| . 

For the PB stochastic process one finds 



-x) 



1 + 
1 + 



q(q - 1) + g(g-l) 2 (g-2) + q(q - l)\q - 2)\q - 3) + 



x 

X 



2\x 2 



3!x 3 



(6) 



7i + 

for the normalized factorial moments [55, and 



K q = q \ x 



l-q 



(7) 



for the normalized cumulants. Here, L^\a) is the associated Laguerre polynomial []60 



The PB stochastic model has been applied to describe the multiplicity distribution and its 
moments in the entire (pseudo)rapidity range [|58] and in its bins [[)9[ in pp collisions at cm. 



energy ranging from 11.5 to 900 GeV. Considering high order moments, it was shown |59[] that 
the results of this approach are close to the NB predictions revealling the stochastic nature of 
particle production and, in particular, of the NB model. Further analysis |35|| , which includes 
intermittency study extending from e + e~ to nuclear collisions^, have shown that the data is 
well reproduced by this model. Note that the predictions of this model are systematically below 
the NB calculations. 

The last form of the multiplicity distribution, considered in this paper, is the recently intro- 
duced generalized negative binomial distribution called the HNB distribution due to the 
type of special H-function used to derive it H6l, 62|. This distribution represents an extension of 
the NB regularity to the Poisson-transformed generalized T-distribution by incorporating some 
perturbative QCD characteristics. Varying the shape parameter k > 0, the scale parameter 
A > and the scaling exponent /i ^ 0, one receives special and limiting cases of the HNB 
distribution. The Poisson and the LN distributions are special cases of the HNB distribution 
in the limit of k —>■ oo with \i — 1 and \i — > 0, respectively. The HNB regularity converges to 
the NB distribution at /x = 1. 

Applied to high-energy data [|6TH64|1, the HNB distribution has been found to agree with the 
data, depending on the parameter \i values being positive or negative^] > 1) or approaching 
zero, for different types of reactions. The HNB regularity with [l > 1 and k = 1 (the Weibull 
law) has been found to describe the data successfully in inelastic pp and pp reactions up to 
ISR energies and in deep-inelastic e + p scattering at HERA energies in the entire rapidity phase 
space |62| as well as in its restricted bins ]63|]. The multiplicity distribution from UA5 data of 
non-diffractive pp collisions at y/s = 900 GeV has been shown to be fitted reasonably well by 
the LN (/i — * 0, k —>■ oo) HNB limit in the full pseudorapidity window and in its symmetric 
bins. 

In e + e~ annihilations, it was found that HNB describes the data below the top PETRA 
energies for \i < — 1 and k = 1 [BE], whereas at high energies the HNB description favour 

IBTl |62|, Bll . The latter, as well as the above mentioned success of the 



/i and large k 



62 



6 See footnote 2. 

7 For /i<0 the following reparametrization of the HNB density has been suggested ]6l| , |6^ ]: (fc, A, «-» (p, a, a) 
with p = l/Vk, a = p/fi and a = In A. With these parameters, e.g. one gets the LN distribution when p = 0. 
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LN limit HNB distribution to reproduce the non- diffract ive UA5 data, is ascribed to the LN 
regularity of the multiplicity distribution^ obtained earlier [^7], [3(| [45| and recently explained 
by a renormalization group approach p2"] . For the LEP1 data /i- value lying between —1.2 and 
—0.6 and k ~ 20 -r- 130 have been obtained. Note that while the error range of the parameter 
\i is found to be small, so that fi is above —2.2 and below zero, the errors for k allow it to vary 
between O(10) and O(10 4 ) fl61 |. The situation does not change when the energy increases up 
to higher than the Z° peak or when the multiplicity distribution in rapidity bins instead that 
in full phase phase space is considered [j6"2| . It is interesting that for central rapidity bins, /i is 
obtained to be positive, < jj, < 1, increasing with enlarging the bin size, while k decreases and 



is of O(10). Recent analysis |64[] of the full phase space of uds-quark jet of the OPAL data |)7 



has shown the error bars for the fc-parameter to lie between 20 and 420, with the central value 
at 54. The parameter \i has been found to be about —0.5. 



One gets |6T| for the HNB-defined normalized factorial moments, 



r(fc + q/fM) 

r(k) 



r(k) 



r(* + i/Ai) 



(8) 



approaching asymptotically for large q to the T-function of the rescaled rank, qj \x. 



3 Comparison with OPAL measurements and discussion 

In Figs. [I] and |^ we show, respectively, the normalized factorial moments and the normal- 
ized factorial cumulants, measured by OPAL || in e + e~ annihilation and compared to a few 
parametrizations (lines) and to the MC |68| (shaded areas). The moments are represented 
in one-, two- and three dimensions of the phase space of rapidity, transverse momentum and 
azimuthal angle calculated with respect to the sphericity axis. 

For the NB law we used the second-order factorial moments or cumulants to compute k and 
then the factorial moments and cumulants of order q > 3, according to Eqs. (jl]) and (@). The 
resulting quantities are shown by the dashed lines. 

From Fig. [I] one can see that in general the NB regularity underestimates the measured 
factorial moments. The deviation is more pronounced in one and two-dimensions for number of 
bins M > 5 (in one projection) and q > 3. Note that these our conclusions coincide with those 
from the analogous investigations |34| of earlier LEP1 results on intermittency [|33|. The better 



agreement between the parametrization and the data we find for low-order (q = 3) moments 
or for those in three dimensions, the cases when the NB predictions are within the MC results. 
Nevertheless, even in three dimensions, the NB values are below the data points at large M 
(small bin size) and high orders. 

The situation becomes more clear when one addresses to the cumulants, namely the genuine 
correlations contributing to the fluctuations, Fig. |[ From one- and two-dimensional cumulants 
one can see that not only for high g's, but even in the case of q = 3 the correlations given by 
the NB model are weaker than those from the data. The same is observed for p?, yxpr and 
Qxpf projections (not shown). The NB cumulants lie far away from the measured ones in 
comparison with the MC predictions being much nearer to the data. Moreover, it is seen that 

8 It is worth to mention the remarkable different shapes of the multiplicity distribution in these two cases, 
namely, a heavy-tailled form with a narrow peak in pp collisions vs. a bell-like shape in e + e" annihilation. 
Nevertheless, recently it was found that the multiplicity distribution in both types of interactions show the 
same behaviour expected from the log-KNO scaling pq| . 
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the discrepancy between the data and the NB results does not begin at some intermediate M, 
as in the case of the fluctuations, but is visible even at smaller M- values (larger bin sizes). 
This observation agrees with the inadequacy of the NB regularity to fit the full phase space 
multiplicity distribution, i.e. at M — 1. The NB parametrization describes the data reasonably 
well in three dimensions, although some deviations are seen for M > 125. 

Using Eq. (fj) we have estimated the parameter A; as a function of M at different q. The 
parameter is found to decrease with increasing M and to depend weakly on q. The values of k 
lie between ~ 0.2 2.8 and ~ 5-^8 at two extreme M values and vary at fixed M with a change 
in q and with the dimensionality of subspaces. The higher is the dimension of the subspace, the 
smaller is the lower bound of the /c-range. These lower bound values are almost independent 
of q. Conversely, the values of upper bound on k show their g-dependence. They are about 8 
at q = 2 and about 5 at q = 3 and 4 regardless of the subspace dimension. According to the 



expectations plU , the observed values of k and its g-dependence do not seem to be related to 
the truncation effect but rather to small cascades. On the other hand, taking into account a 
long enough cascade at the e + e~ collision energies considered here, this conclusion confirms that 
NB encounters difficulties in a reasonable and consistent description of the measurements. It is 
interesting to note that the values of k obtained are close to those found |54| in the MNB-type 
analysis of the multiplicity distributions in restricted rapidity intervals in e + e~ annihilations at 
the Z° peak g§. 

Contrary to the NB, the LN regularity overestimates the data regardless of the dimension- 
ality or the type of the variable used. The dotted line in Figs. [I] shows how the LN predictions 
compares the measurements. The cumulants are calculated from the factorial moments us- 
ing their interrelations 0. The smaller the bin size is, the larger the difference is. The LN 
parametrization describes the data quite well for order q = 3 only. 

These findings are in agreement with the earlier studies of LN fits to factorial moments 
in e + e~ annihilation |34|. The LN distribution overestimates high multiplicities^, which leads 
to an overestimate of the fluctuations and genuine correlations as shown here. In contrast to 
the studies of the full- multiplicity parametrization p^,RU|, the deviations are found for all bin 



sizes and not only for intermediate ones. This is in contrast with the above mentioned multi- 
jet (perturbative) effects and indicates significant contribution of the non-perturbative stage 
dynamics, i.e. soft hadronization, to the formation of fluctuations and correlations. This agrees 
with recent theoretical studies . Violation of the Gaussian law of Eq. (|3|) have already been 
observed in nuclear [37,40, 4^] and hadronic [6£| interactions. 

The predictions of another, the MNB regularity, are shown in Figs. [I] and ^| by the solid 
lines. The cumulants are calculated using Eqs. (f|) and (||), while the factorial moments are 
derived from their relationships with cumulants ||. The parameters r and A have been fixed 
at the values r = 0.91 and A = —0.71, the best values found to describe at least the third order 
cumulants. The only parameter depending on the bin size is k extracted from K 2 . The value of 
the A parameter is near to that found in multiplicity studies [p8|,^7l,^9|-|5T|, while the parameter 
r has the value r > |A|, in contrast to the r < |A| inequality obtained from the analysis of the 
full multiplicity distributions. From the Figures, one can conclude that, in general, with the 
parameters obtained, the MNB regularity describes the data well, although underestimates the 
latter at small bins and the parameters |A| and r have an inverse hierarchy. 



9 See, however, Q. 

10 A better description of the multiplicity distribution in the case of smaller number of jets in e + e~ annihilation, 
namely for two-jet events, compared to the inclusive sample, has been also obtained by OPAL using the NB 
regularity p7fl. 
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The obtained exceeding of the value of r over the A-value is expected if one applies the 
MNB regularity not to negativ e-particle distributions but to all-charged-particle onesQ see e.g. 
DELPHI publication PS|. In this case, two effect are contributing: the number of sources 



increases with increasing the width of the phase space bin, while a neutral cluster decay gives 
0, one or two particles hit in the given bin ||57|p|. The latter effect is taken into account by 
implementing an additional parameter [54]. This is not the case when corrected formulae (|5|) 
are used, therefore the change of the inequality between r and |A| has a physical meaning to 
be investigated. 

The resulted factorial moments and cumulants of the multiplicity regularity given by the 
PB process, are represented in Figs. [I] and § by the dashed-dotted lines. The parameter x is 
extracted from the second order cumulants K 2 and defines all the higher-order moments and 
cumulants, given by Eqs. @ and (JF]), respectively. For both quantities the PB predictions 
are seen to lie lower than those from the above considered parametrizations. In all the cases 
the PB calculations underestimate the data. The difference between the PB predictions and 
the data is in contrast to earlier lower-energy e + e~ parametrization of the factorial moments 
in the rapidity subspace |B5|, where the PB process has been found to explain the data while 
underestimating higher-order moments. 

No curves predicted by the HNB regularity are shown in the Figures, since we would like 
to estimate the regions of the parameters /i and k of this approach. Combined analysis of the 
factorial moments and cumulants, the latter being calculated as combinations of the factorial 
moments 0, show that, assuming k to be positive, the /i-parameter is obtained to be either 
positive or negative. The factorial moments and cumulants are found to be sign-changing 
functions of k for negative /i, in particular at — 1 < fi < 0, so that one can find more than one 
fc-region which satisfies the measurement sPl In this case we take the largest /c-value, above 
which no sign-changing behaviour is seen and the calculation fits most of the data. For fi > 
the factorial moments and the cumulants decrease with increasing k. 

At fixed k, the absolute value of the parameter /i decreases with the number M of bins. The 
parameter is found to be limited in the interval 0.3 < < 1.8 when one considers the quantities 
under study of order q = 2,3. For q > 3, one faces problems reaching the highest measured 
values of the moments and cumulants and, in the case of the highest q = 5 cumulants, fitting 
their lowest (negative) values. This narrows the interval of /i down to 0.6 < < 1.6. The 
values of the shape parameter k depends on /i, but are found to be ^ 30. The smallest k is 
obtained to be 1.2 for /j, < and 0.1 for /x > 0. The larger the value of |/x| is, the smaller is the 
interval of k. For example, if at |/i| > 1.3 the values of k are only of a few units, ~ 2 4- 5, then 
for 0.5 < \fi\ < 1.0 these values lie between 2 and 20. The parameter k slightly depends on 
the order q, increasing with q for negative /i and decreasing for positive ones. This parameter 
depends also on the number of bins, being a decreasing function of M. Note that these two 
properties are similar to those observed here for the NB regularity. 

Comparing the results obtained to the HNB studies |JT],^,|65| of the LEP data and, par- 
ticularly, to that of the uds-quark jets |64| (see also Sect. ||), one can see similarities as well as 

11 See, however, footnote 5. 

12 It is interesting to note that being limited to like-charged particles, a study of particle bunching is less 
dependent on correlations induced by charge conservation (and partly by resonance production) , in addition to 
the above-mentioned advantage of such a study to be less affected by the energy-momentum constraints fl57f . 

13 Our observation contradicts the property of the factorial moments and cumulants shown (6l],|6^] to oscillate 
around zero at > 1 and not e.g. at /i > — 1. This disagreement can be assigned to different regions of the 
parameter k, found to be large (k —> oo) in the case of the full-multiplicity distribution studies |]6l| , [65| while 
having finite values in our investigation (vide infra). 
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important differences. Similar to those studies we find that the parameter /x can be positive as 
well as negative and does not exceed 2 in its absolute value. Moreover, this parameter shows the 
same behaviour with bin size, and k tends to obtain large values |62,E3j]. The main difference 



from the HNB studies is that the region \i « is excluded by our investigation, so that the 
LN law is not a suitable one. This conclusion agrees with the OPAL observation |27] and the 
above discrepancy between the data and the LN predictions. It is also worth to notice that 
(i) negative /i's can be also used to describe the multiplicity distributions in restricted bins, 
in addition to the positive ones found recently ||64j| , and (ii) the values of the /c-parameter are 
less than 30 and do not tend to infinity. All this shows that to describe the hadroproduction 
process correctly, one needs a more complicated scenario to be realized than those leading to 
the regularities discussed here, even generalized to the HNB case. One could not e.g., find /i to 
be only in the interval < /i < 1 as it would expected due to our findings for the LN regularity 
(/x ~ 0) to overestimate the measured fluctuations and correlations while the NB law (// — > 1) 
underestimates them. 



4 Summary and Conclusions 

To summarize, variuos regularities of the multiplicity distribution of charged particles are stud- 
ied in restricted phase-space bins of e + e~ annihilation into hadrons at the Z° peak. The study 
is based on recent high-statistics results || on multidimensional local fluctuations and gen- 
uine correlations obtained with multihadronic OPAL sample by means of normalized factorial 
moments and cumulants. Such parametrizations as the log-normal and negative binomial dis- 
tributions, modified and generalized versions of the negative binomial law, and the generalized 
stochastic birth process with immigration are considered. For the first time these parametriza- 
tions, being most common in the field of multiparticle production, are examined with genuine 
high-order correlations. 

All the parametrizations are found to give a reasonably good description of low-order fluctu- 
ations while they do show deviations for the high-order fluctuations and correlations, especially 
at small resolutions. Some discrepancy can arise from a bin-size dependence of the measured 
multiplicity distribution and from its truncation. However, in our consideration, the influence 
of these effects is minimized since the models are based on the measured second-order factorial 
moments and cumulants, which carry most of the information given by the multiplicity dis- 
tribution. Moreover, the effects mentioned arise mostly according to low statistics that is not 
the case for the data considered here, even at higher orders. To note is also that even with 
low statistics data, a simultaneous analysis of multiplicity distributions in different bin sizes 



and the corresponding factorial moments, carried out in e + e annihilation [34] and in nuclear 



reactions |40|1 , does not show any sensitive influence of the finite statistics to the results. 



From the study presented, one concludes that genuine high-order (q > 3) correlations have 
to be taken into account when the hadroproduction process is modelling, in particular by 
a multiplicative law for particle distributions. Indeed, since all the parametrizations used are 
essentially based on the average multiplicity and the two-particle correlations, the discrepancies 
between the predictions and the measurements indicate multipaiticle character of bunching of 
hadrons. This could be considered as a reason why all the regularities give a good description 
of fluctuations and correlations at order q = 3. Our conclusion confirms the OPAL result 
from Ref. ||, on which we are based here and which shows the important contributions of 
many-particle correlations to the dynamical fluctuations by the decomposition of the factorial 
moments into lower-order cumulants. This our finding is also in agreement with the observation 
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of DELPHI |7| shown the measured angular cumulants not to be reproduced by small genuine 
correlations given [p ,[l5| by perturbative QCD. 



A self-similar nature of multihadron production is another issue of the above study. All the 
regularities result from the particle-production process of the stochastic nature. Therefore their 
capability to show the intermittent behaviour seen in the data can be attributed to their branch- 
ing self-similar nature. The better agreement between the regularities and the measurements 
found in three dimensions, where QCD cascading is expected to be fully developed, stresses 
the essential self-similarity of the particle production mechanism. However, the discrepancies 
obtained show that a suitable hadroproduction model seems to be more sophisticated than 
those giving the parametrizations considered in the present paper. 
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Figure 1: Factorial moments of order q = 3 to 5 as a function of M 1//D , where M is the number 
of bins of the D -dimensional subspaces of the phase space of rapidity, azimuthal angle, and 
transverse momentum, in comparison with the predictions of various multiplicity parametriza- 
tions and two Monte Carlo models. The data and the Monte Carlo predictions are taken from 
Ref ®. 
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Figure 2: Cumulants of order q = 3 to 5 as a function of M X I D , where M is the number of bins 
of the D -dimensional subspaces of the phase space of rapidity, azimuthal angle, and transverse 
momentum, in comparison with the predictions of various multiplicity parametrizations and two 
Monte Carlo models. The data and the Monte Carlo predictions are taken from Ref. f^j. 
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